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1 Introduction 

General relativity and quantum mechanics together imply that space-time structure at the 
Planck scale is described by noncommutative geometry [1]. There have been various at- 
tempts to study gravity theories within the noncommutative framework [2, 3]. This has led 
to a Hopf algebraic description of noncommutative black holes [4, 5] and FRW cosmologies 
[6]. A large class of such black hole solutions, including the noncommutative BTZ [7, 8] 
and Kerr black holes, exhibits an universal feature where the Hopf algebra is described 
by a noncommutative cylinder [9], which belongs to the general class of the k- Minkowski 
algebras [10-13]. For the purpose of this paper, we shall take the noncommutative cylinder 
and the associated algebra as a prototype for noncommutative black holes. 

The study of quantum field theories in the background of black holes has led to the dis- 
covery of interesting features associated with the underlying geometry, such as the Hawking 
radiation and black hole entropy. In the noncommutative case, the black hole geometry is 
replaced with the algebra defined by the noncommutative cylinder. In order to probe the 
features of a noncommutative black hole, it is useful to analyze the behaviour of a quantum 
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field coupled to the noncommutative cylinder algebra. Scalar field theories have been ex- 
tensively studied on K-Minkowski spaces [14-17, 19], which has led to twisted statistics and 
deformed oscillator algebra for the quantum field [16-18]. Theories on the noncommutative 
cylinder lead to quantization of the time operator [20, 21] 

The quantum field theories defined on noncommutative spaces are highly nonlocal and 
in order to gain further insight into their behaviour, it is essential to simulate their be- 
haviour through numerical analysis. To this end, it is necessary to approximate the infinite 
dimensional noncommutative cylinder algebra with a suitable truncated finite dimensional 
matrix algebra, belonging to the general class of fuzzy spaces [22-24]. 

Field theories on noncommutative geometries are inherently as mentioned earlier non- 
local leading to mixing of infrared and ultraviolet scales. This, in turn, is responsible 
for new ground states with spatially varying condensates. Many non-perturbative studies 
have established that noncommutative spaces, such as the Groenewold-Moyal plane and 
fuzzy spheres, allow for the formation of stable non-uniform condensates as ground states. 
Exploring the implications of the non-local nature of field theories is very important in 
many areas of quantum physics. [25, 26] 

Since different phases are intimately connected with spontaneous symmetry breaking 
(SSB), the role of symmetries in noncommutative geometries themselves is subtle. This 
issue is important in 2D because the Coleman-Mermin- Wagner (CMW) theorem states 
that there can be no SSB of continuous symmetry on 2-dimensional commutative spaces. 
There is no obvious generalization of the CMW theorem for non-commutative spaces, since 
the theorem relies strongly on the locality of interactions. Noncommutative spaces admit 
non-uniform solutions (in the mean field) and one can ask the question what happens to 
the stability of these configurations. Non-uniform condensates naturally have a infra-red 
cut-off for the fluctuations. This cut-off softens the otherwise divergent contributions of 
the Goldstone modes [27-34]. 

This paper is organized as follows. In Sec. 2 we give a brief introduction to fuzzy black 
holes motivated by an earlier analysis of the noncommutative deformation of Banados- 
Teitelboim-Zanelli (BTZ) black holes. In Sec. 3 we set up the algebra describing the 
noncommutative cylinder which is suitable for numerical simulations. In Sees. 4 and 
5, we provide the action for the scalar field and the spectrum of the Laplacian on the 
noncommutative cylinder respectively. Sec. 6 exhibits the phase structure and the novel 
stripe phases which are generic to fuzzy spaces. We exhibit the crucial differences of these 
different phases and analyze our results in Sec. 7. 



2 Fuzzy black holes 

We briefly summarize the essential features of a noncommutative black hole which is useful 
for our analysis. In the commutative case, a non-extremal BTZ black hole is described in 
terms of the coordinates (r, (/>, t) and is given by the metric [7, 8] 

/ 2 7-2 \ / 2 t2 \ -1 / / \ 2 

ds'=^M-^-^yt'+^-M + ^ + —^^ dr' + r'[d4>-—,dtj , (2.1) 
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where < r < oo , — oo<t<oo, < < 2tt , M and J are respectively the mass and 
spin of the black hole, and A = — is the cosmological constant. In the non-extremal 
case, the two distinct horizons r± are given by 

An alternative way to obtain the geometry of the BTZ black hole is to quotient the 
manifold AdS^ or SL(2, R) by a discrete subgroup of its isometry. The noncommutative 
BTZ black hole is then obtained by a deformation of AdS-s or SL{2, R) which respects the 
quotienting [4]. In the noncommutative theory, the coordinates r, <j) and t are replaced by 
the corresponding operators f, (j) and Z respectively, that satisfy the algebra 

[i, e'^] = ae'^ [f, t\ = [r, e'^] = , (2.3) 

where the constant a is proportional to ^^/(r^ — r^). We shall henceforth refer to (2.3) as 
the noncommutative cylinder algebra. Furthermore, t denoting the operator corresponding 
the the axis of the cylinder, it will be therefore identified as the operator Z also in the 
following sections. 

It may be noted that the operator r is in the center of the algebra (2.3). In addition, 
it can be shown easily that g"^'^**/'^ belongs to the center of (2.3) as well. Hence, in any 
irreducible representation of (2.3), the element g"^'^**/" is proportional to the identity, 

e"2^**/° = e^^n, (2.4) 

where 7 € mod (27r). Eqn. (2.4) implies that in any irreducible representation of (2.3), 
the spectrum of the time operator t, or Z, is quantized [9, 20, 21] and is given by 

spec t = na — — , n € Z. (2.5) 

In what follows we shall set 7 = without loss of generality. 

The noncommutative cylinder algebra (2.3) belongs to a special class of the K-Minkowski 
algebra and it appears in the description of noncommutative Kerr black holes [5] and FRW 
cosmologies [6]. We shall henceforth consider (2.3) as a prototype of the noncommutative 
black hole. 




3 Noncommutative Cylinder Algebra 

The NC cylinder is defined by the relation 

Z, e'* 



ae 



(3.1) 



where Z is hermitian and e**^ is unitary. As mentioned earlier, the operator Z corresponds 
to the axis of the cylinder, and therefore to the time operator t of the black hole. Since we 
are interested in simulations, we have to discretize the above NC cylinder. 
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In the rest of this paper, we will work with a = 1 without loss of generality, since the 
simple scaling Z — )• Z/a can scale a away in the commutation relation (3.1). 

For this purpose consider the spin J irreducible representation (IRR) of the SU (2) Lie 
algebra, given by 

= 2X3, [X±,X3] = tX±. (3.2) 

Since the operator Z generates rotations around the axis of the cylinder, it can be identified 
with X3. But when we use the finite dimensional representations of SU{2) we cannot 
implement (3.1) with unitarity for e*"^. For this purpose, we decompose X+ as product of 
a unitary and a Hermitian operator as given by 

X+ = e''^ R (3.3) 

In (3.3), e^"^ is unitary, and i? is a positive Hermitian, necessarily singular, matrix which 
commutes with Z (and is thus diagonal). Using (3.2) and the fact that R commutes with 
Z, we have 

[Z, X+\ = [Z, e^<^] R = e'^ R (3.4) 

Since R is singular, it can not be inverted. However a partial inverse R can be found such 
that RR = P, the projector such that 1 — P projects on the kernel of R. Thus we get 

[Z,e''^] P= e'^ P (3.5) 

To find a representation for R and e^'^ , we can look at 

i?2 = = P- L^iLs + 1) (3.6) 

which commutes with Z = L3. Remember that in the usual representation of angular 
momentum, L3IZ, m >= m\l, m > with |?7i| < I. Shifting the indices from to 2Z + 1 = 2 J, 
we have j = m + Z + 1, leading to 

X_X+\l,m > = + -m(m + l)]|/,m >= [{l + l/2f - {m + l/2f]\l,m > 
X_X+|J,j > = [J^-{J-jf]\J,j >=j{2J-j)\J,j > 

There is only one hermitian positive solution to this equation which takes the form 



Rij = ,Ji{2J - i) 5,,j (3.7) 

which is diagonal as expected, and whose null space is along the top state |J, 2J >. As a 
result, P = 1 - I J, 2 J >< J, 2J|, and 

R=Y^ [i{2J - i)]-i/2 I ^ >^ ^1 

4 = 1 

where the sum stops at i = 2 J — 1 so that there is a zero in the last position on the diagonal. 
It is now possible to deduce the first 2 J — 1 lines of the unitary matrix e*"^ from (3.3): 



X+l J,i > = Vj(2J~j)\J,j + 1 > 



X+\J,j > = e^^ R\J,j >= ^j(2J-i)e'^ \J,J > 



\JJ >= + l >, j < 2J. 
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Eq. (3.3) yields no equation for the last column which is instead determined from its 
unitarity. The columns 1, • • • , 2J — 1 of e'*^, given by | J, j + 1 >, form an orthonormal set, 
as expected for a unitary matrix. Then the last column will be a vector orthogonal to all 
these vectors and thus can only be proportional to | J, 1 >. After normalization that still 
leaves a C/(l) freedom so that 

= 6,,,+i + e'^ 6,,i6,,2j (3.8) 

where /? can be any real number. For (3 = 0, e^'^ is just a circular permutation of length 
2J. 

4 The Action on the fuzzy cyhnder 

We will first construct an action for a hermitian scalar field Define TrO = Tr{POP). 
This trace Tr is equivalent to integrating over the whole cylinder in the continuum limit. 
We also need the derivatives and dz- They are: 

= [z,$] (4.1) 

dz^ = e-^'^[e^'^,$] (4.2) 

Then, apart from J-dependent normalization factors, a naive form of the action can be 
chosen as: 

S = fr(\ |2 + I e-''^[e'^,<^] ^ + T/($)) (4.3) 

where V{^) is the potential which can be taken to be of the form, 

for a hermitian field 

This action has a problem of instability. The source of this comes from Tr($^) = 
Tr{{P^)^^{^P)) which cannot contain any quartic (nor cubic) term for the variable $2J2J- 
This makes the theory unstable with respect to this variable. The simplest cure is to insist 
that this term is not a degree of freedom of the theory and constrain it to zero. To keep 
the set of fields an algebra, we choose to also set to zero the last row ^2Ji and column ^i2j 
of the field. As a result the hermitian field <I> now only has (2 J - 1)^ de grees of freedom, 
and $ = P^P. 

4.1 Dimensional reduction 

With this new choice of the field the action becomes 

S = Tr(^P \ [Z,P<^P] + P I e-''^[e'*,P^P] p P+ V{^)^ (4.5) 

= Tr (\ [PZP,$] p + I [Pe"f'P,^] p + y($)) (4.6) 

which can be rewritten simply as the action for a hermitian matrix in a (2 J — 1) x (2 J — 1) 
matrix algebra of reduced dimension: 

S = Trj>(^\ p +1 [?^,^.] |2+ l/($)^ (4.7) 
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where J' = J — 1/2 is the reduced angular momentum, while e**^ and Z are the matrices 
obtained from e**^ and Z by removing the last line and column. For e*''^, it is equivalent 
to setting e!"^ — ?> in its 2J'-dimensional expression (3.8). As for Z, it is therefore the 
2J' X 2J' diagonal matrix obtained from Z by removing its top eigenvalue J — 1/2: 

Z = Diag(- J + 1/2, - J + 3/2, • • • , J - 3/2) = Diag(- J', - j' + 1, • • • , j' - 104.8) 
^ Zij = {-J' -l + i)6i^j. (4.9) 

Note that Z and Z are defined by their commutation relation (3.1) and only appear in the 
action through 5^ as a commutator. As a result, they are only defined up to a translation 
by a matrix proportional to the unit, and thus Z = Diag(l, • • • , 2 J') is another possible 
choice. 

Although e^'t' is not unitary, the equation 

< m! \[Z ^ e'^'t'Wm > = {rn! — m)6m',m+i = Sm',m+i =< m'\e^'f'\m > if m < 2J' 
= =<m'\[Z,e^<t>]\2J' > if m = 2 j' 

shows that Z and e*"^ do satisfy the commutation relation (3.1) for a = 1 

4.2 Expressing the Kinetic term 

Using the reduced action (4.7), the kinetic term K(^) then takes the general form 

_ 2J' 

K{<^)=Tr(\[Z,<^] |2 + \[e^<P,^] = ^ [(^ _ j-)2|$^ .p + |ci>._^ . _ ci>.^.^^j2] ^ (4.10) 

where we have introduced new entries $oj = ^i2j'+i = set to zero to simplify the 
expressions. 

The kinetic term can be further reordered as 

2J' i-1 2J 

^(^) = Yl - + 4 - S^,N - 5,,i]|$i,f + ^(2 - 5i,N - S^,i)\^ii\^ - (4.11) 

i=l j=l i=l 
2J'-1 i-1 2J'-1 

1=1 j=i 1=1 

Note in this expression that the last two sums are not over all possible indices, but omitting 
the highest one. Furthermore, in the last two terms, is coupled to both <l>j+ij+i and 

^^-lJ-l 

The expression of the potential (4.4) is already known, being the same expression as 
for the fuzzy sphere (see e.g. [27]). 

4.3 The action on a fuzzy cylinder of radius r 

The cylinder is also parametrized by its radius r. According to (2.3), r commutes with 
both Z and e**^. It can therefore be considered as a pure number in the non-commutative 
cylinder algebra. 
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The radius will appear as a simple scaling in the action. The volume of the cylinder 
Tr{l) depends linearly on r, so the action should have an overall scale of r. The derivative 
along the axis dz does not scale with r, whereas the angular derivative scales like 1/r. 

As a result, the action on a fuzzy cylinder of radius r is given by: 

S = rfr (^\ [Z,$] p + I e-'^[e^'^,$] p + J (4.12) 



r Trj, ( 4| |2 +1 [ei<t>,^ p + ) (4.13) 



^ I IZ.0I r + I \e,^<P.<P\ r + 1/(0! 
5 Spectrum of the Laplacian on the fuzzy cyhnder 

The Laplacian comes from the kinetic term (4.10) of the action. After integrating by parts^, 
we get 

K{^) = rr(-[Z,$][Z,^>] - [?^\$][e^,^>]) = rr($[Z, [Z, ^>]] + [e^,$]]) 
so that naively 

= [z, [z, ^>]] + [?^\ [?^, $]] = /:|^> + 

where £2, resp.>C^, resp. C_, is the adjoint action of Z^ resp. e*"^, resp. e**^ . Note however 
that this Laplacian is not hermitian due to e**^ not being actually unitary in the second 
term. To make it hermitian, the latter term must be symmetrize. The Laplacian now is: 

£2$ = [Z, [Z, <5]] + ^([e^\ [e^, $]] + [7^, $]]) 

= Cl^ + ]^{C_C^^ + C^C_^) (5.1) 
= + + ^[[e^,e^^], + + ^[Diag(-l,0,--- ,0,1),$] 

The Laplacian derived above in (5.1) is for a cylinder of radius one. The action (4.13) 
on a cylinder of radius r shows how the Laplacian scales with the radius r, yielding 

£2$ = J-£|$ + i(£ £ $ + £ £ $). (5.2) 
r 2 

Since we will be interested in the entropy of the free field, let us now look for the 
eigenvalues of the Laplacian. 



^ Since commutators work as derivations and the trace as integration, there is a direct equivalent to 
integration by part given as 

rr($[L, *]) = Tr([L, - Tr([L, $]*) = -Tr{[L, 
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5.1 Symmetries of the Laplacian: 

Because the Laplacian is expected to be invariant with respect to Z-axis rotations, it must 
commute with Lz- This is quite obvious on the expression of the Laplacian since: 

- e^'i' has axial momentum +1 and will therefore raise total axial momentum by one 
whether multiplied on the left or on the right 

- conversely e^'^ has axial momentum —1 and will therefore lower total axial momen- 
tum by one whether multiplied on the left or on the right 

and therefore overall the axial angular momentum is conserved by C_C^ or C_C^. Taking 
into account that we have hermitian eigenmatrices, we can deduce that they must be a mix 
of +m and —m axial momentum matrices. This means that we can take the ansatz 

2J'-m. 

$ = $^((f) + ^i,{d), with ^m{d) = ^ |i >< i + m|, < m < 2/ - 1 (5.3) 

i=l 

for the eigenmatrices, where the vector d = ((ii)i<i<2j'-m is the unknown to be determined 
by the eigenmatrix equations. 

Using this ansatz (5.3) and the hermiticity of $ and the Laplacian, the eigenmatrix 
equation for the Laplacian reads simply 

C^<^m{d) = X'^mid) = $m(Ad) ^ M,^d = Xd 

where the matrix Mm translates the (linear) action of the Laplacian on <^m((i) to d: 
jC'^^mid) = ^miMfnd). This 2J' — m square matrix is now evaluated. 

Note that the Laplacian has another less obvious symmetry, it is invariant under the 
replacement \i >^ \2J' — i >. We will not need to use it in the following calculation, but 
the eigenvector equations in the following should be (and we checked that they actually 
are) invariant under that symmetry. 

5.2 The matrix Mm. 

Piece of Mm coming from C"^: In this case, since by construction, $m is an eigenma- 
trix of jCz with axial angular momentum m, 

Cl^m{d) = m'^<^rn{d) = <^m{m^d)- 
Piece of Mm coming from £|^>C j: 

(J) = di/:{^C_}{\i >< i + m\) 

i 

= di {{x^,x^}\i >< i + m\ + \i >< i + m|{x+, x_} — 

i 

>< i + m\x- — x^\i >< i + m\x-^.) 

= '^^i ((1 ~ ~ '5i,2J'/2 + 1 - <5i+m,l/2 - Si+m,2J'/'^)\i >< i + m\- 

i 

(1 - '^j,2J')(l - '5j+m,2J')N + l><i + m+l\- 
(1 - <5i,i)(l - 5i+m,i)\i -l><i + m-l\) 



-8- 



where x- = e^'t' lowers the index of the ket, whereas = x^__ raises it. At this point in 
the calculation, it is necessary to distinguish the cases m = 0, for which all the Kronecker 
can be 0, and m ^ where only half of them can. 
For m = 0, we find 

C{^C_}<Pm.{d) = ^di ((2 - - 5i^2j')\i >< «| - (1 - 5i,2j')N + 1 >< ^ + 1|- 

i 

{I - 6i,i)\i - 1 >< i - 1\) 

/ 1 -1 \ 

-1 2 -1 (0) 



(0) 



v 



2 -1 
-1 1 



d) 



For m ^ 0, on the other hand, 
C^^C_}^m{d) = Y,di{{2- <5i,i/2 - Si+m,2.r /2)\i ><i + 



m\ 



(1 - 6i+rn,2j')\i + 1 >< i + m + 1| - (1 - 5i^i)\i - 1 >< i + m - 1|) 
/ 3/2 -10 \ 



-1 2 -1 



(0) 



(0) 



V 



d) 



-.2-1 

-1 3/2 J 

Expression of Mm Putting together the results from the last two paragraph, and using 
the Laplacian (5.2) on a cylinder of radius r, we get 



/ 1 -1 
-1 2 -1 



Mn 



(0) 



(0) 



V 



. 2 -1 

-1 1 



/ 3/2 -1 
-1 2 -1 



"1 



\ 



(0) 



(0) 



V 



■. 2 -1 
-1 3/2 J 



, m/0 (5.4) 



These matrices are similar to the ones obtained for the Laplacian on a one-dimensional 
lattice and can actually be diagonalized without much difficulty, taking good care to re- 
member that Mm is a matrix of dimension 2J' — m. 

5.3 Spectrum of Mm 

The way to evaluate the spectrum is to write explicitly the eigenvector equation. Then, let 
d = (di) be an eigenvector for the eigenvalue A, the eigenvector equation takes the form 

dn+i = (2 - A)(i„ - dn~i, 2 < n < 2J' - m (5.5) 
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plus boundary equations at each end n = 1,2J' — m, which are different for m = and 
m 7^ 0. 

The sequence defined by this linear induction formula with constant coefficients can be 
determined by looking at its characteristic equation^ q'^ — {2 — X)q + 1 = 0. Reparametrizing 
the eigenvalues as 

A = 2 - 2cos(6') = 4sin2(6'/2), (5.6) 

this equation has the simple solutions q = exp(ibi^). Therefore, the sequence has the 
general form 

dn = 7^(/3e*"^) (5.7) 

where /3 is a complex constant to be determined by the first two terms of the sequence. 

To simplify the algebra, it is convenient to extrapolate do and d2j-m+i and rewrite 
the boundary equations between them, di, and d2j-m- Furthermore, since the eigenvector 
is defined up to an overall constant, let us choose do = 1. 

Now we must look at the two cases separately. 

Spectrum of Mq In this case, the boundary equations can be seen to take the form 
Similarly, at the other boundary, 

d2j'-i = d2j'. (5.8) 

In particular, do = di = 1. Plugging those initial values in the general form (5.7) of 
the sequence, we get 

and therefore 

c^n = ^;^cos((n-l/2)^). (5.9) 
To get the eigenvalues, it only remains to enforce the last constraint equation (5.8): 
cos((2J'- 1/2)61) = cos((2J' + 1/2)6*) ^ sin(2 J'6l) sin(6l/2) = 0^9 = kir/2J', 0<k< 2J', 
which, according to (5.6) gives eigenvalues: 

Ag = 4sin2(A;7r/4J'), < < 2J'. (5.10) 
A corresponding eigenvector is then given in (5.9) up to an overall constant as: 

d= (cos((?i - 1/2)W2J')),<„<2J' ■ (5-11) 



^The linear space of sequences satisfying the induction formula (5.5) is of dimension 2, parametrized by 
the two initial values of the sequence. The idea is to find a basis of this linear space in the form of two 
geometric sequences of the form (g"). q must then satisfy a (quadratic) characteristic equation, and the 
sequence we want to express is a linear combination of these two geometric sequences. 
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Spectrum of M^, m ^ In this case, the boundary equations can be seen to take the 

di = 2dQ 



form 

f d2 = (3/2 - X)di 



[ da = (2 - X)di - do 
And similarly, at the other boundary, 

2(i2J'+l-Tn = d2J'-m- (5.12) 

In particular, do = 1 and di = 2. Plugging those initial values in the general form (5.7) of 
the sequence, we get 

/ 7^(/3) =1 , .cos(0)-2 



7^(/3e*^) = 2 sm{e) ' 



and therefore 



sin(6l) ^ ^ sin(6l) sin(6l) ^ ^ 

To get the eigenvalues, it only remains to enforce the last constraint equation in (5.12). 
Denoting N = 2J' — m, and using the trigonometric relation 

sin(a) -sin(6) = 2cos((a + 6)/2) sin((a - 6)/2), 

we get: 

2(sin((A^ + 1)61) - sin(iV6')) = sin(iV6l) - sin((iV - 1)9) 
44> 2 cos((A^ + 1/2)61) = cos((iV- 1/2)61) <^ 3tan(iV6l) tan(6l/2) = 1 
<^ tan(7V6') = tan((7r - 6')/2)/3. (5.14) 

which can be seen graphically to have N solutions, one in each interval [kir/N; (A; + l)7r/A^], 
< k < N . These solutions must be determined numerically though. The eigenvalues are 
then deduced from Eq. (5.6). The corresponding eigenvector is then given by (5.13). 

For large matrices N ^ 1, it is possible to find approximate solutions since tan(x) can 
be well approximated near and 7r/2. This yields: 

• For 6* < vr, or A; < N, tan((7r - 6l)/2) ~ 2/6* » 1. Therefore, NO = kn + 7r/2 - pk 
with ^ 1. The equation then becomes: 

1 2N 37r , , 

^" 3(i^vr + vr/2-p,) ^^^"2iv(^ + ^/^^- 

• For 61 ~ TT, or /c ~ N, tan((7r - 6')/2) ~ (vr - 6l)/2 < 1, and therefore, N9 = kn + pk 
with Pk <^ I. The equation then becomes: 



A^vr — kn — Pk N — k 

m ^^'^"^^ 

which is a small number, as expected, since /c ~ A. 



- 11 - 



6 Numerical Simulations and Results 



6.1 The numerical scheme 

The model defined by the action Eqs. (4. 13,4.4) which we want to simulate has three param- 
eters (//, c, r) plus the matrix size J. The goal is to explore the parameter space for various 
phases of The simulations are carried out using the "pseudo-heat bath" Monte-Carlo 
(MC) algorithm [31, 32] to reduce the auto-correlation along the MC history. 

The field should also be allowed to explore the whole phase space and not remain 
trapped in local minima. To this end, an over-relaxation method, first suggested in [30], 
is also used. Let us introduce 5$($ij) the dependence of the action on the field entry (^ij 
when the field takes the value It is a fourth degree polynomial. Therefore the equation 
S^{^ij) = S^{^ij = a), which has an obvious solution = a, can be factorized into a 
degree three polynomial which always admits at least one real solution. The overrelaxation 
method consists in replacing the field entry = a by one of these real solutions, thereby 
moving the field in a different region of the phase space. 

A crosscheck is also used to verify that the field probability distribution of our Monte- 
Carlo runs are consistent. Let us split the terms in the action according to their scalings 

S{(l)) = 52(-/.) + 54(0) with SiiXcP) = X'Si{<p). 

Then one can define a modified partition function 

Z{X) = [ [d0]e-^(^'^) = / [dcP]e-^'^^(^)-^'^^('i>) (6.1) 



A-^ / [d^]e-^(^), V = A0, (6.2) 



where is the number of degrees of freedom in the field (p which appear in the integration. 
Evaluating 

d\n{Z) 



dX 



= -2 < 52 > -4 < ^4 > from (6.1) 

A=l 

= -N from (6.2) 



yields the check originally due to Denjoe O'Connor [35]. 

< 52 > +2 < 54 >= N/2 (6.3) 
In all simulations, this identity (6.3) is always satisfied to better than 1% relative error. 
6.2 The phase structure 

The temperature (T) is regulated by varying the parameter fj,. 

• fi <^ 1 corresponds to low temperatures when the fluctuations are small. In this case, 
the minimum of 5 gives the most probable configuration of the phase. In Eq. (4.13), 
it is possible to minimize the action by minimizing separately the kinetic term, so 
that $ oc 1, and the potential term so that <I> = \J —^ijlc 1, and this phase is therefore 
known as the uniform phase. 
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• At high temperatures, fi ^ 1, the thermal fluctuations lead the system to the disorder 
phase <I> ~ 0. 

• At intermediate temperatures, the competition between the action and the fluctu- 
ations give rise to new phases called the non-uniform or stripe phases. These new 
phases are speciflc to non-commutative spaces. Various numerical studies have con- 
flrmed the existence of these phases [27-34] on the fuzzy sphere. A non-commutative 
cylinder will also exhibit the non-uniform phases. However, due to the non-trivial 
topology of the cylinder (the flrst homotopy group being non-trivial), one can have 
a more complex phase structure described below. 

For example there can be stripes going around the cylinder, or parallel to its axis. 
These two phases can be distinguished by their overlap with the operators Z, e^'^, and e^'^^ 
respectively. Stripes going around the cylinder will have non-zero overlap with the operator 
Z. While a configuration of stripes along the axis will have overlap with e^'^ and e*"'^^. We 
present our results in the following subsection. 

6.3 Example numerical runs 

For a given choice of = 7, c = 0.36, and r = 1 the simulations are done for various 
values of /i. The various phases discussed above can be characterized by the observables 
m„ = Tr($), niz = Tr(<I>Z), rux = Tr{^e^'^). A finite m„ with (mz,mx) ~ characterizes 
the uniform phase. On the other hand, (m„, nix) ~ with non-zero rriz characterizes stripes 
going around the cylinder. Stripes along the cylinder characterized by (jn^^ m^) ~ with 
non-zero m^. 

For /i = —35.1, the data of a run are shown on Figure 1, and, as expected, we observe 
the uniform phase. 

For fj, = —20.0, we observed the phase with stripes going around the cylinder. This 
is verified on the histogram of the observed values of mu,mz plotted in Figure 2. It is 
clear from the figure that the average value of niz is finite while the average value of niu is 
vanishingly small. 

Figure 3 shows the system in the disorder phase where niu, mz,mx all fluctuate around 
zero. 

We did not observe the phase with stripes going along the cylinder as a ground state 
for any choice of ^ for r ~ 1. One can expect to observe this state for very small r when 
the second term of the kinetic term in the action (4.13), which suppresses this state, is 
made subdominant. For a very small radius r = 0.01, c = 36. and fi = —3.6 x 10'^, this 
phase appears as meta-stable in Figure 4. This phase is stable w.r.t small fluctuations. 
Only large fluctuations, which occur less frequently, can destroy such a state. 
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Figure 1. rriu, niz, ra^ vs Monte Carlo history ioi N = 7 , ji = —35.1, c = 0.36, and r = 1. 
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Figure 2. Histogram of H{mu) and H{mz) N = 7, fi = —20.0, c = 0.36, and r = 1. 
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Figure 3. Histogram of H{mu) and H{mz) N ^ 7, c ~ 0.36, and r = 1. 
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Figure 4. rriu, mz,mx vs Monte Carlo history for = 7, /i = —3.6 x 10'^, c = 36., r = 0.01. 



7 Conclusions 

In this paper, we have considered a finite dimensional representation of the noncommutative 
cylinder algebra, which makes it fuzzy. We study scalar field theory in the background of 
this algebra both analytically and using numerical simulations. 

The action of the scalar field on the fuzzy cylinder contains a kinetic as well as a 
potential term. The kinetic term leads to the Laplacian (5.2) on the fuzzy cylinder. We 
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have analyzed the symmetries of the Laplacian and have obtained an algebraic equation 
(5.10,5.14) describing the corresponding spectrum. 

In the numerical simulations of the scalar field with a generic potential we find, as 
expected in noncommutative theories, novel stripe phases breaking rotational symmetry. 
But they have some differences with the usual stripes on Moyal spacetimes. These are also 
stable due to topological features arising in this fuzzy geometry. A similar stability has 
been seen in the 0(3) model on fuzzy spheres [36]. 

The fuzzy cylinder algebra considered in this paper is valid for non-extremal BTZ black 
holes, without any further assumptions. In addition, different forms of the BTZ metric 
related by coordinate transformations are equivalent classically, which at the algebraic 
level are expected to be related by automorphisms. 

While the noncommutative cylinder algebra first arose in the context of the BTZ black 
hole, subsequently it has been shown to be of more general relevance, appearing in diverse 
backgrounds such as for Kerr black hole [5] and FRW cosmologies [6]. In addition, the 
near-horizon geometry of a large class of black holes contains an AdS^ factor [37] . 

Thus, a large class of noncommutative black holes are described by a noncommutative 
cylinder algebra. The fuzzy cylinder algebra derived from it can therefore be used to define 
a fuzzy black hole. From general considerations [1], we know that such black holes can 
arise at the Planck scale. Our results provide a first glimpse about the phase structure of 
a quantum scalar field theory in the background of a fuzzy black hole at the Planck scale. 
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